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CN ■ Abstract We discuss the dynamical effects of bulk viscosity and particle creation on the early 

^ . evolution of the Friedmann -Robertson -Walker model in the framework of open thermody- 

^!^ namical systems. We consider bulk viscosity and particle creation as separate irreversible 

^ ■ processes. Exact solutions of the Einstein field equations are obtained by using the "gamma- 

^^ ■ law" equation of state p = (7 — l)p, where the adiabatic parameter 7 varies with scale factor 

of the metric. We consider the cosmological model to study the early phases of the evolution 
of the universe as it goes from an inflationary phase to a radiation -dominated era in the 
presence of bulk viscosity and particle creation. Analytical solutions are obtained for par- 



bJO tide number density and entropy for all models. It is found that, by choosing appropriate 

functions for particle creation rate and bulk viscous coefficient, the models exhibit singular 
Cf^ ■ and non-singular beginnings. 
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1. Introduction 

The study of particle creation in the relativistic cosmological models has drawn the atten- 
tion of a number of authors. The first theoretical approach of the particle creation problem 
were investigated by Prigogine et al.^'^. They showed that the second law of thermodynamics 
may be modified to accommodate flow of energy from the gravitational field to the matter 
field, resulting in the creation of material particles. This leads to the reinterpretation of 
the stress tensor in general relativity, which now involves a time-asymmetric term depend- 
ing on the rate of creation of particles. In other words, the process of particle creation out 
of gravitational energy is basically an irreversible phenomenon, capable of explaining the 
entropy burst in the expanding universe. A detailed study of the thermodynamics of the 
particle creation with changing specific entropy have been discussed by Lima et al.^, Calvao 
et al.^, Lima and Germano ^, Lima et al.^, and Zimdahl et al.^. Johri and Desikan ^ studied 
Friedmann-Robertson- Walker (FRW) models with particle creation. 

After the discovery of the accelerating universe the particle creation was reconsidered to 
explain it and got unexpected results. The particle creation pressure, which is negative, might 
play the role of dark energy component. Lima and Alcaniz ^, and Alcaniz and Lima ^^ dis- 
cussed FRW models with matter creation and analyzed the results through the observations. 
It was shown that the models with particle creation are consistent with the observational 
data. Zimdahl et al.^"^ tested the particle creation with SNe la data and got the result of 
the accelerating universe. Yuan Qiang et al.^^ studied the models with adiabatic particle 
creation and showed that the model is consistent with SNe la data. 

On a phenomenological level particle creation has been described in the literature in terms 
of a bulk viscous stress. Padmanbhan and Chitre ^^ discussed the role of the bulk viscosity 
in the entropy production in an expanding universe. Another peculiar characteristic of bulk 
viscosity is that it acts like a negative energy field in an expanding universe ^^. The basic idea 
was that the bulk viscosity (particle creation) contributes at the level of the Einstein field 
equations as a negative pressure term. Barrow ^^'^^ introduced this idea in the framework 
of new inflationary scenario. However, Prigogine et al.^'^ pointed out that the bulk viscosity 
and particle creation are not only two independent processes but, in general, lead to different 
histories of the evolution of the universe. Sudharsan and Johri ^'^ , Brevik and Stokkan ^*, 
and Triginer and pavon ^^ discussed the effect of bulk viscosity on the cosmological evolution 
of open thermodynamic systems which allow for simultaneous particle creation and entropy 
production. Their investigations reveal that the production of specific entropy is independent 
of the nature of the creation rate and depends only upon the coefficient of the bulk viscosity. 
This suggests that the bulk viscosity and particle creation are not only independent pro- 
cesses, but in general they lead to different histories of cosmic evolution. Desikan '^^, Singh 
and Beesham ^^, Johri and Pandey 22, Singh et al.^^, and Singh and Kale ^^ have studied 
the role of particle creation and bulk viscosity in isotropic and anisotropic models. 

It has been suggested in the literature that dissipative processes in the early stage of the 
expansion of the universe may well account for the present high degree of isotropy. A detail 
study of string -driven inflationary universe in terms of effective bulk viscosity coefficient was 
studied by Barrow ^^. Gr0n ^^ and Maartens ^^ presented exhaustive review on cosmological 
models with non-causal and causal thermodynamics respectively. Several authors ^7-29 g^g_ 



gested that the bulk viscosity can drive the universe into a period of exponential expansion 
(inflation). This is really the case, as the effect of bulk viscosity in an expanding universe is 
to decrease the pressure making the total pressure negative. Chimento et al.^" found that a 
mixture of a minimally coupled self-interacting scalar field and a perfect fluid is unable to 
drive the accelerated expansion, while a mixture of dissipative fluid with bulk viscosity and a 
minimally coupled self-interacting scalar field can successfully drive an accelerated expansion. 

In the standard model, the history of the universe begins with the radiation phase and 
then evolves to the present matter -dominated era. In order to overcome some of the difficul- 
ties met by standard model, Guth ^^ proposed an inflationary phase and this would happen 
prior to the radiation -dominated phase. In general, the field equations are solved separately 
for the different epochs. However, Some authors have tried to solve the field equations in a 
unified manner. Madsen and Ellis ^^ presented the evolution of the universe for inflation- 
ary, radiation and matter-dominated phases in a unified manner by assuming gamma (7) 
of "gamma-law" equation of state p = (7 — l)/9 as a function of scale factor of the FRW 
metric. Later on, Israelit and Rosen ^^'^^ used a different equation of state to describe the 
transition from pre-matter to radiation and then radiation to matter-dominated phase in a 
unified manner. 

In a similar way, Carvalho ^^ studied flat Friedmann-Robertson- Walker (FRW) model in 
general relativity by using the "gamma-law" equation of state where 7 varies with cosmic 
time to describe the early phases (inflation and radiation) of the evolution of the universe 
in a unified manner. Therefore, it it not realistic to assume 7 as a constant throughout 
the history of the universe. We can obtain a reasonably realistic model if we assume the 
universe evolves through the epoches each of which 7 is constant. Singh et al.'^^, Singh '^^ 
studied flat viscous FRW model and showed that the viscous fluid might drive the present 
acceleration of the universe. Recently, Singh and Beesham ^^ have studied FRW model with 
the particle creation for the early phases of the evolution of the universe and discussed the 
kinematic tests. 

The aim of this paper is to extend Carvalho's work ^^ to include the theory of particle 
creation and viscosity. We discuss the dynamical effects of bulk viscosity and particle creation 
on the early evolution of the Friedmann -Robertson -Walker model in the framework of open 
thermodynamical systems. We consider cosmological model with bulk viscosity and Particle 
creation as separate irreversible processes. Exact solutions of the Einstein fleld equations 
are obtained by using the "gamma-law" equation of state p = {"f — l)p, where the adiabatic 
parameter 7 varies with scale factor of the metric. We consider the cosmological model to 
study the evolution of the universe as it goes from an inflationary phase to a radiation - 
dominated era in the presence of bulk viscosity and particle creation. Analytical solutions 
are obtained for particle number density and entropy for all models. It is found that, by 
choosing appropriate functions for particle creation and bulk viscous coefficient, the models 
exhibit singular and non-singular beginnings. 

2. Basic equations 

We start with the homogeneous and isotropic fiat Freidmann-Robertson- Walker (FRW) 



line element in the units SvrG = 1 and c = 1 

ds'^ = dt^ - R^{t)[dr'^ + r^{de^ + sin^ed^^)], (1) 

where r, 0, and (j) are dimensionless comoving coordinates and R is the scale factor. The 
Einstein's field equations are given by 

where Tij is the effective energy momentum tensor of the cosmic fiuid in the presence of 
creation of particle and bulk viscosity, which includes the creation pressure term, pc and the 
bulk viscous stress, 11 and may be defined as "^^ 

Tij = ip + P + Pc + Il)uiUj - {p + Pc + ^)gij (3) 

where p is the energy density, p is thermodynamical pressure and Ui is the four velocity vector 
satisfying the relation u^Ui = — 1. The creation pressure, pc is associated with the creation 
of particle out of the gravitational field ^'^. The bulk viscous pressure, 11 represents only 
a small correction to the thermodynamical pressure, it is a reasonable assumption that the 
inclusion of viscous term in the energy -momentum tensor does not change fundamentally 
the dynamics of the cosmic evolution. 

In context of 'open' system with adiabatic creation, the non-trivial Einstein's field equa- 
tions for a fiuid endowed with matter creation and viscosity can be written as 

3;^ = P^ (4) 

R R , ^^N / s 

2;^ + ;^ = -(p + Pe + n), (5) 

where an overhead dot denotes derivative with respect to cosmic time t. For adiabatic particle 
creation, the pressure pc assumes the following form ^'5:20 

{p^p)V dN 

where V = R^ \s the 3-space volume and N is the particle number. In models with adiabatic 
creation, the balance equation for the particle number density n, where n = N/V is given by 

4,20 

N _n R _ T{t) 

N~n^^R~~^' ^^' 

where T{t) denotes a source term, i.e., the particle creation rate which will be positive (F > 0) 
or negative (F < 0) depending on whether there is production or annihilation of particles. 
In this work we take the simple phenomenological expression of particle creation rate ^ 

F(t) = S^nH, (8) 



where the parameter /? is defined on the interval [0, 1], which is assumed to be constant and 

H = R/R = {1/3){V/V) is the Hubble parameter. 

Using (7) and (8) into (6), the particle creation pressure reduces to 

p, = -l3{p + p). (9) 

Equations (4) and (5) lead to the continuity equation 

p + 3{p + p)H = -3{pc + U)H. (10) 

The conventional bulk viscous effect in a FRW universe can be modelled within the 
framework of non-equilibrium thermodynamics proposed by Israel and Stewart^'^. In this 
theory, the transport equation for the bulk viscous pressure 11 takes the form 



n + rll = -3CH - ^ 



n T c 



:n) 



where the positive definite quantity C stands for the coefficient of the bulk viscosity, T is the 
temperature of the fluid, and r is the relaxation time associated with the dissipative effect. 
We consider here the role of bulk viscosity in framework of equilibrium thermodynamics. 
Although equilibrium thermodynamics has some shortcomings, in contrast to the extended 
irreversible thermodynamics, the treatment of bulk viscosity under equilibrium thermody- 
namics in the present paper is justified for the following reasons. 

Provided the factor in the square bracket is small, one can approximate the equation (11) 
as a simple form 

n + ril = -3CH, (12) 

Assuming r = C/p as taken by Maartens ^^, and noting that the energy density p is propor- 
tional to the square of the scalar expansion 9, we have 11 = — C^[l + H/^^]. It is seen that we 
can take the equilibrium value for 11, namely, 

n = -3CH, (13) 

provided 11/^^ <^ 1. In the particle creation models, 9 happens to be a rapidly increasing 
function of time, and hence the condition 11/^'^ <^ 1 holds. The entropy equation is assumed 
to be of the form 

where S is the entropy and 9 = 3{R/R) is the expansion scalar. In terms of specific entropy 
per particle, o", Eq.(14) can be written as 

We see that the production of specific entropy per particle is independent of the nature of 
r(t) and depends only on the nature of bulk viscosity. If (^ = 0, we have a =constant. 



3. Solution of field equations 

In order to solve the field equations, we suppose that the pressure p and energy density 
p are related through the "gamma-law" equation of state 

P={l- 1)P- (16) 

In general, the value of 7 is taken to be constant and lying in the interval < 7 < 2. 
But our aim in this paper is to let the parameter 7 depends on scale factor R to describe 
the early phases, inflationary and radiation -dominated evolution of the universe in a unified 
manner. We assume, following Carvalho ^^ that, the functional form of 7 as 

A A{R/R,f + {a/2){R/R,Y 
^^ ' 3 A{R/R,Y + {R/R,Y ' ^ ' 

where A is a constant and 'a' is free parameter related to the power of the cosmic time t 
during the inflationary phase. Here, R.^ is a certain reference value of R. The function 7(-R) 
is defined in such a manner that when the scale factor R is less than i?*, i.e., when R « R^,, 
an inflationary phase (7 < 2a/3) can be obtained and for R >> i?* we have a radiation- 
dominated phase (7 = 4/3). The expression of 7(-R) in Eq.(17) is an increasing function of 
R. In the limit i? — )■ 0, 7(i?) = -^. Thus, 1 is the maximum value of 'a' for an inflation epoch 
to exist. As 'a' approaches to zero we have an exponential inflation (7 = 0). Therefore, a 
must lie in the interval < a < 1. 

Using (9), (13) and (16) into (10), we flnd 

p + 3j{l-(3)pH = 9CH\ (18) 

which can be written as 

where a prime denotes derivative with respect to the scale factor R. 
Equation (4) can be rewritten as 

p' 2H' , , 

7 = IT- P°> 

Using (4) and (20) into (19), we finally get 

if' + |(l-/3h(H)f = 5i. (21) 

Equation (21), involving H and C, admits solution for H only if (^ is specified. In the ho- 
mogeneous models, C depends only on time and therefore we may consider it as a function 
of the universe energy density. According to literature developments [see refs., Belinskii and 
Khalatnikov ^^, Barrow ^^ and Maartens ^^], we assume that the bulk viscosity coefficient 
depends on p via a power -law of the form 

C = Cop"^, (22) 



where Co is a non -negative constant and m{> 0) is numerical constant to be specified later. 
It is standard to assume the above law in the absence of better alternatives. In what follow 
we solve Eq.(21) by taking the various physical assumptions on m 

3.1 Solution with m = 



In this case, we get the constant coefficient of bulk viscosity, i.e., C = Co a-i^d therefore, Eq. 
(21) now reduces to 



^'.f(i-,«h(/.)f^i. 



Substituting (17) into (23) and solving, we get 



^ ^3Co f [A{R/R,r + {R/R,r](^-P) 
Co + — j^ dR 



(23) 



(24) 



where Cq is a constant of integration. We now solve Eq.(24) for two early phases of the 
universe viz. inflationary and radiation -dominated phases, respectively. 
For inflationary phase {R « i?*), Eq(24) gives 



H 



Co 



+ 



3Co 



{R/R,)^(^~'^) 2a(l-/3)' 
and for radiation-dominated phase {R » R^), we have 

Co , 3Co 



H 



+ 



A(i-/3)(i?/i?,)2(i-/3) 4(1-/3)' 



(25) 



(26) 



From above we find that the Hubble parameter is constant when the viscous term dominates 

in both the phases and therefore the solutions are de -Sitter type where the scale factor varies 

as i? oc exp(3Co/2a(l — /?)) or i? oc exp(3Co/4(l — /3)) and density is finite. The universe has no 

singularity. When the particle creation dominates we get the singular model as i? oc t^'"-(^^P) 
ori?ociV2(i-/3). 

When Cq ^ 0, after some mathematical manipulations between the constants , Eq. (25) 
gives 

'exp(3fit)-(l-/3)' 



R* 



2aB 



3Co 



where B is another constant. The Hubble parameter has the form 



H 



3Co 



l-(l-/3)exp(-^i 



-1 



2a(l - /3) 
The energy density and particle pressure are respectively given by 



27Co^ 



4a2(l-/3)2 



l-(l-/3)exp(-^t 



(27) 



(28) 



(29) 



Pc 



9/3Co' 



2a(l-/3)2 
The particle number density is given by 



l-(l-/3)exp(-%t 



1 -2 



R 

n = noi\ -^ 



-3(1-/3) 



(30) 



(31) 



We note that the effect of particle creation is measured by the parameter /?. The particle 
number is given by 






3/3 



(32) 



In the above expressions the subscript "Oi" refers to the present observed values of the 
parameters during inflationary phase. It is observed that the universe starts from a non- 
singular state. We observe that for p > 0, (3 must be in the interval < /3 < 1, which 
is also the condition for the expansion of the universe. It shows thatp starts with a finite 
value at t = and ends up with a finite value at t — )■ oo. It is interesting to note that the 
inflationary solution can be obtained with or without particle creation and this is due to 
constant bulk viscous coefficient. In the absence of both bulk viscosity and particle creation 
one get the power-law inflation, i.e., R oc t^'"', which exhibits singular model. We also find 
that N increases with R increases. For (3 = 0, N would remain constant throughout the 
evolution of the universe. 

In order to get an expression for a we need an additional equation of state involving T. 
For the sake of simplicity we take the ideal gas equation p = nkf,T with k^, the Boltzmann 
constant, which we use just as an approximation since ideal gas lack bulk viscosity. Inserting 
this value into Eq.(15) for 7 = 2a/3 we get 



a 



9Coh 
(2a - 3) 



t + l, 



(33) 



where I is an integration constant. 

On the other hand, when the universe is dominated by radiation, i.e., when R >> R^,, 
Eq.(26) gives 

^ 4i?i fexp(%)-(l-/3)' 



^X 2(1-/3) 



yl(i-/3) [ 3Co 

where Bi is a constant of integration. The Hubble parameter has the form 

3Co 



H 



l-(l-/3)exp(-^t 



-1 



4(1-/3) 
The energy density and particle pressure are respectively given by 



P 



27Co^ 



16(1-/3)2 



l-(l-/3)exp(-^t 



(34) 



(35) 



(36) 



9/3C, 



'l-(l-/3)expf-^t'' 



Pc = -4^^^[l-(l-/3)exp(^-^tJJ . (37) 

The particle number density and correspondingly particle number are given by Eq.(31) and 
(32) where the scale factor is defined by Eq.(34). The specific entropy per particle is found 
to be 

a = 9Cokbt + li, (38) 

where li is constant of integration. 

In this phase the physical interpretation is similar to the case of the inflationary phase. 
The bulk viscosity avoids the singularity. For the expansion of the universe, we must have 
0</3< 1. 

The deceleration parameter, which is defined as q = —{RR/B?), varies from q = [a(l — 
/?) exp(— 3CoV2) — 1] for infiationary phase lo q = [2(1 — /3) exp(— 3^0^/2) — 1] in radiation 
-dominated phase. We find that q = [a(l — /3) — 1] in inflationary phase and ? = (1 — 2/3) in 
radiation phase as t = 0. However, it is g = —1 in both phases as t — t- oo, which shows the 
inflation during the late times of evolution of the universe. 

Now, we study the model in the limit a — t- and in this case Eq. (24) becomes 

H [MRIR.f + 1] "-"• = Co + ^ / \MRIR.f + lf^-'\^^ ,3^) 

Li it 



Again, in the limit of very small R {R « R^,), we get 

R = R^, exp 



exp(2fit)-Co 



3Co/2 



(40) 



which means that the universe has 'superinflationary' expansion due to the presence of viscous 
term. The effect of particle creation is negligible. It has flnite dimension as t — t- — oo. The 
Hubble parameter has the value H = exp{3Cot/2) which shows that the expansion is driven 
by viscosity. 
The energy density is given by 

p = 3exp(3Cot) (41) 

The energy density is flnite at t = and hence there is no physical singularity. These solutions 
help to solve several cosmological problems like flatness, horizon, monopole, etc. associated 
with standard model. The particle creation pressure is zero. The superinflation is due to 
either by the bulk viscous coefflcient or vacuum energy density {p = —p). Thus the cosmo- 
logical constant may also be considered to give rise the exponential inflation in the absence 
of creation pressure. We also observe that the particle density decreases whereas the number 
of particle increases exponentially during this phase. 

In the limit of very large R, when the universe enters to radiation-dominated phase, 
we obtain the same physical expressions for different physical parameters as in the case of 
< a < 1. We see that the deceleration parameter varies from g = — lati? = Otog= (1 — 2/3) 
for radiation-dominated phase. 



3.2 Solution with m / 

Using (4) and (22), Eq.(21) now reduces to 

H'^l(r-0Un)!^.'^S^. (42) 

Substituting (17) into (42) and integrating, we have the solution, for all m ^ 1/2, which is 
given by 



H('^m^^)[A{R/R,)^ + (i?/i?,)"](i-/3)(2™-i) 
(2m - l)3-"+iCo f dR 

' 2 y i?[A(i?//?,)2 + (i?/i?,)«](i-/3)(2-"i)' ^ ""^ 

where Ci is a constant of integration and < a < 1. We discuss the solution for two early 
phases of the universe viz. inflationary and radiation-dominated phases respectively. 
For inflationary phase {R « i?*), the expression for Hubble parameter is given by 

jji2m-l) ^ 1 /^4^ 

[Ci(i?/i?*)"(2™-i)(i-/3) + (3('"+i)Co/2a(l-/3))]' ^ ' 

The energy density, creation pressure and bulk viscous coefficient in terms of scale factor for 
this phase are respectively given by 

p = 3[Ci(i?/i?*)"(2'»-i)(i-/3) + (3™+iCo/2a(l - /3))]2/(i-2'»)^ (45) 

Pc = -2a/3[Ci(i?/i2,)"(2'"-i)(i-^) + (3™+iCo/2a(l - /3))]2/(i-2m)^ (45) 

( = 3"'Co[Ci{R/R,f^^"'-^^^^-'^^ + (3™+iCo/2a(l - /3))]2™/(i-2m). (47) 

An important observational quantity is the deceleration parameter q. In the case of infla- 
tionary phase, the deceleration parameter can be expressed as a function of the scale factor 
as 

^ (a(l - /3) - l)Ci(i?/i?.)"(2"^-^(i-/^) - (3'»+iCo/2a(l - /3)) . . 

^ [Ci(i?/i?,)'^(2™-i){i-/3) + (3™+iCo/2a(l -/?))] ^ ^ 

The sign of q indicates whether the model inflates or not. The negative value of q {q < 0) 
describes the acceleration of the universe whereas the positive one (q > 0) describes the 
decelerating universe. In the absence of viscous term, it becomes q = a(l — /3) — 1. This 
shows that the universe decelerates for a(l — /3) > 1 and accelerates for a(l — /3) < 1. Since 
< a < 1, in the absence of creation pressure, i.e., /3 = and Ci > 0, we get the negative 
value of q showing the inflation in the model. 
For radiation -dominated phase (R » i?*), we have 

j^(2m-l) ^ 1 (^g^ 

10 



In this phase the energy density, particle pressure and bulk viscous coefficient in terms of 
scale factor are given by 

p = 3[Ci(i2/i?*)2(2™-i)(i"'^) + (3™+iCo/4(l - /3))]2/(i-2™), (50) 

p, = -4/3[Ci(i2/ii*)2(2™-l)(l-'3) ^ (3^+1^^/4(1 _ ^))]2/(l-2m)^ (5^) 

C = 3™Co[Ci(i?/i2*)^(^""^)(^-^) + (3™+^Co/4(l - /3))]2™/(i-2™). (52) 

The deceleration parameter, in this case, in terms of scale factor can be written as 

^(2m-i)(i-/3)(i _ 2/3)Ci(i?/i2*)2(2™-i)(i-/3) - (3™+^Co/4(l - /3)) 



[^(2"^-l)(l-/3)Ci(i?/i?,)2(2m-l)(l-/3) + (3m+l(o/4(l _ /J))] 



(53) 



In the absence of viscous term, we have q = {1 — 2/3). This shows that the model decelerates 
for < /? < 1/2 and accelerates for 1/2 < /3 < 1 in the absence of viscous term. If Ci > 0, 
we observe that q is positive for (i?/i?*)2(2™-i)(i-/3) > (3'"+iCo/4(l - /3)A(2m-i)(i-/3)(^^)^ 
q is negative for (i?/i?,)2(2™-i)(i-/3) < (3™+iCo/4(l - p)A^^"'-^^^^-P^Ci) and g = for 
(it:/i?,)2(2'»-i)(i-/3) = (3™+iCo/4(l - /3)^(2m-i)(i-/3(^^)^ Again, if Ci = 0, we get q = -1, 
which shows the acceleration of the universe. 

We discuss the more general solution of the above models for the following three cases 
depending on the value of m for inflationary and radiation -dominated phases. 

3.2.1 Solution vi^ith m > 1/2 

In this case, we observe that Eqs. (44) and (49) give constant value of Hubble parameter 
as i? — )■ 0. The model is similar to those of Murphy '^'^ , which has initially the viscosity 
-dominated steady state behavior. The bulk viscosity appears to be the effective mechanism 
to remove the initial singularity. Let us discuss this case by taking m = 1 as the bulk viscosity 
is proportional to the energy density. In this case, Eqs. (44) and (49) for inflationary and 
radiation -dominated phases, respectively can be explicitly obtained as 

^' ^ ^ ^ "^'^^^ + . 5° .. lnR = t + to. (54) 



a(l-/3) \RJ 2a(l-/3) 

Ci fRY^'-^^ 9Co 



2il-P)\Rj +4(1-3^^"^ = * + *- ^''^ 

where to is the integration constant and can be adjusted to zero. In general, it is not pos- 
sible to express R explicitly as a function of time. For sufficiently small R, the second 
term in equations (54) and (55) dominate over the first, containing the viscous term. Thus, 
we get exponential expansion of the form R = exp(i7*t), where -ff* = 2a(l — /3)/9Co or 
H^ = 4(1 — /3)/9(o, which represents the singularity free model and all solutions approach 
to the de Sitter state with the expansion rate H = H^. Such a non -singular behavior is 
exhibited only in the presence of bulk viscosity. Such type of solution has been obtained 
by Murphy (1973), where he attributed the viscosity effect to gravitons production in the 

11 



graviton -gravitons scattering. Since H^ '^ Co ' ^° ^^^ smaller is the value of Co the bigger is 
the inflation rate. During inflation both p and C are finite. The deceleration parameter has 
the value q = —1. The above solutions (54) and (55) can also be directly obtained from the 
general solution (44) or (49) after some manipulation. 

On the other hand, if the first term dominates we get the power -law expansion R oc 
ii/"(i-/3)^ which gives the singular model and the expansion is driven by particle creation. 
We observe that the effect of viscous coefficient becomes negligible. The model shows sin- 
gularity and the expansion continuously slows down but never reverses. Thus, we see that 
this model evolves from a de Sitter state as t — t- — oo to zero curvature Friedmann state as 
t — 7- +00. At t = 0, the scale factor and energy density are finite, where as the scale factor 
tends to zero and the energy density becomes infinite as t — )■ —00. The viscosity removes the 
initial singularity at finite past i.e. moves it to the infinite past. The energy density and bulk 
viscosity decrease as time passes. Similarly, for radiation -dominated phase where the bulk 
viscosity dominates we get R oc ex.p{Ht), where H = 4(1 — /3)/9Co- In the absence of viscous 
term, we have R oc t^/^li-/^)^ p oc t~^ and C oc t~^. The deceleration parameter has the value 
q = 1, which shows the deceleration of the universe. 

3.2.2 Solution with m < 1/2 

In this case, Eqs.(44) and (49) for inflationary and radiation phases can respectively be 
rewritten as 



H 



H 



Ci(i2/i?*)-''(i-2m){i-/3) ^ (3{'"+i)^Q/2a(l - /3)) '-'™ . (56) 



Ci(ii/i?.)-2(i-2™)(i-^) + (3(™+i)Co/4(l - /3)) '"''" . (57) 



In this case, the universe evolves into a viscosity dominated steady state era. If the coefficient 
of bulk viscosity decays sufficiently slowly, the late epochs of the universe will be viscosity 
dominated, and the universe will enter a final infiationary era with steady state characteris- 
tic. This case is distinct from Murphy's case or other values of m with 2m > 1. This can 
also be explained from (44) and (49). If the coefficient of bulk viscosity decays slowly, i.e. 
if 3™+^Co/2a >> Ci{R/R^Y^'^"^~^\ there is exponential expansion R oc exp{Hot). Thus, at 
any finite proper time in the past, the curvature is finite. The viscosity has removed the 
initial singularity. On the other hand if Ci{R/ R^,)"'^'^"''^^' >> 3™'^^Co/2a, we get power -law 
expansion R oc t^'°-(^^l^) and the effects of viscosity are negligible. Similarly, we can discuss 
the behavior of solution in radiation -dominated phase. 

3.2.3 Solution with m = 1/2 

The solutions obtained in subcases (3.2.1) and (3.2.2) are not valid for m = 1/2. Let us 
consider the case m = 1/2 and in this case C = CoP^ ■ Now, Eq.(42) becomes 



H' + 



3(1 - PMR) 
a 
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- = 0. (58) 



where a = (3v3Co/2) • Substituting (17) into (58) and integrating, we have the solution for 
Hubble parameter 

[AiR/R,)^ + {R/R,)-]i^-^)- ^ ' 

where C2 is the constant of integration, li H = H^ for i? = i?*, a relation between constants 
is given by 



Equation (59) further integrated to give 

[A{R/R,f + (i?/fi,)"](i"/^) 

Jla+l 



dR = C2t (61) 



The constant of integration has been taken as zero for simplicity. The scale factor for infla- 
tionary and radiation -dominated phases are, respectively given by 

and 

/o2\ (1-/3) 

i?[2(i-/3)-"l = (^^ j [{2(1 - /3) - a}C2t] , (63) 

The corresponding Hubble parameters are given by 

and 

^' [2(l-ffl-.| i f"^* 

The above solutions show the power-law expansion of the universe for a(l — /3) > a or 
2(1 — /3) > a. The energy density, particle creation pressure and coefficient of bulk viscosity 
for inflationary phase are respectively given by 

3 1 

(66) 



^ [a(l-/3)-Q]t2' 

2a^ 1 ._. 

Pc = --rn a\ T72' (^'^ 

[a(l — p) — a\ t-^ 

and 

_ \/3Co 1 , . 

Similarly, the solutions of p, pc and C for radiation dominated phase have the following forms 

^=[2(l-'^)-a]^' ('') 
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4/3 1 .„^. 

^- = - [2(l-/3)-a] ^' ^'°) 

and 

C= ^^° 1 (71) 

^ [2(1 - /3) - q] t ' ^ ^ 

Singular solutions are obtained for n = 1/2 with power-law expansion in both phases. The 
evolution begins from singularity with a Friedmann leading behavior and increases in size 
in the course of time with a decrease of energy density. The bulk viscous coefficient also 
decreases with time. 

From (59), a unified expression of the deceleration parameter for both inflationary and 
radiation phases can be expressed as a function of scale factor as 

_(!_«_ 2I3)A{R/R,)^ + [a(l - /3) - 1 -a]{R/R,r ^^^^ 



[A{R/R,f + {R/R,Y] 

Therefore, q varies from q = [a(l — (3) — a — 1] for R « R^, to q = {1 — a — 2/3) for radi- 
ation phase. The deceleration parameter is positive for [a(l — (3) — a — 1] > 0, negative for 
[a(l — f3) — a — 1] < and g = for a(l — /3) — a = 1 during inflationary phase. Similarly, in 
radiation -dominated phase the universe decelerates for 2(1 — /?) — a — 1>0 and accelerates 
for 2(1 - /3) - a - 1 < 0. 

We now study the solution in limit a — )■ 0. In this case, Eq. (59) becomes 

[A(i?/i?,)2 + l](i-/3)' ^'''> 



which can be written as 






The constant of integration has been taken as zero for simplicity. In the limit of very small 
R, the scale factor is given by 

i?° = -^t-\ (75) 

we observe that a > and C2 > lead to contraction. As t — t- — cx) , we find that i? — t- 0. 
The model starts from infinite past with zero proper volume. Thus, for a = the universe 
is infinitely old and we have inverse power -law. Again, the radiation dominated phase is 
described by the same solutions as obtained for < a < 1. 

From Eq.(73), a unified expression for deceleration parameter can be given in terms of 
scale factor as 

_ ^i_a_2f3)A{R/R,f-{l + a) ^^^^ 



A{R/R,Y + 1 

Therefore, q varies from q = —(1 -|- a) for infiationary phase to q = {1 — 2(3 — a) for R » R 
as expected. 
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4. Conclusion 

We have discussed the role of bulk viscosity in the framework of the particle creation mecha- 
nism in homogeneous and isotropic flat FRW model. Analytical solutions for the scale factor, 
energy density, particle creation pressure, bulk viscous coefficient and number particle density 
have been derived for two early phases of the evolution of the universe with varying equation 
of state parameter 7. We have found that some of the models do not have the initial sin- 
gularity due to the exponential inflation or even superinflation expansions. However, some 
models have the singular solutions. 

The concept of viscosity term has been used in a generalized form. We have analyzed 
the consequence of the inclusion of such a dissipative term in both inflationary and radiation 
-dominated phases. The solution of the models are qualitatively similar in both phases. We 
have observed that the bulk viscosity change the behavior from that of the perfect fluid model 
with particle creation. We have seen that the introduction of viscosity term in to the equa- 
tion of Friedmann cosmology does not exclude automatically the appearance of singularity. 
Within the discussed class of models, singular and non -singular solutions appear. 

It is evident from the three cases discussed in section 3.2 that for 2m > 1 and 2m < 1, 
the strong energy condition p + 3p > 0, where p = p + Pc — 3C-ff , is not satisfied for Ci < 0. 
It means that except for Ci positive the energy conditions are violated throughout the evo- 
lution. But for Ci > and 2m, < 1 the energy condition , which is equivalent to i? < is 
satisfied at the initial stage of expansion and violated at later stages, whereas for Ci > 
and 2m > 1 the energy condition is violated at the beginning and satisfied at later stages 
of evolution. When 2m, = 1, the solution is (59), where there is no de Sitter exact solution 
and the general solution is a power -law form, all energy conditions are satisfied. In this 
case there is a big -bang type singularity, where the evolution begins from singularity with a 
Friedmann leading behavior. In case subsection 3.2.2, the weak energy condition p > and 
dominant energy condition p + p > 0, are always obeyed, but the strong energy condition 
p + 3p > is violated at early times, which allows the avoidance of singularity. Therefore, 
there is no problem with the Hawking -Penrose ^^ energy conditions for A; = and 2m = 1, 
because it is satisfied throughout the evolution. But the energy conditions are satisfied for 
part of this period either at initial stages or at later stages according as 2m, < 1 or 2m > 1. 
We also observe that similar type of energy condition dp)/dp < 1 is violated at later stages 
or initial stages of evolution according as 2m > 1 or 2m, < 1 , but there is no problem with 
2m = 1. 

In the case of constant coefficient of bulk viscosity the solutions to the field equations can 
be expressed in an exact exponential form. The universe starts from a non -singular state, 
characterized by constant and finite initial values oi R , H and p. In case where the coeffi- 
cient of bulk viscosity is proportional to energy density, we get a result which is similar to 
the solution obtained by Murphy. We observe that the bulk viscosity is capable of removing 
the cosmological singularity and the model approach to de Sitter universe. The solution is 
quite different in the case when the C, oc p^'"^ . In this case, the solutions correspond to power 
-law infiation in inflationary phase where as power -law expansion in radiation -dominated 
phase. The energy density and bulk viscosity decrease and tend to zero for large time. The 
universe starts from a singular state. We have observed that Murphy's model is the only 
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special case of these general set of solutions. The effect of bulk viscosity is more prominent 
at the beginning of the universe, where the expansion scalar and energy density are quite 
large. The bulk viscosity are negligible for non-relativistic and ultra-relativistic fluids but 
are important for the intermediate cases. As we can analyze from this work that the effect 
of bulk viscosity is negligible in late time. It affects only during the early universe. In this 
paper the bulk viscosity have been investigated mainly to search the non-singular models 
as motivated by Murphy ^^ and Maartens ^^. The particle creation rate increases the rate 
of expansion in each model and exhibits the singularity. The expansion rate depends on 
the parameter /?. Therefore, it is thought that the negative pressure caused by matter cre- 
ation may play the role of a dark energy and drives the accelerating expansion of the universe. 
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